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Abstract , 

^Results  of  an  investigation  into  nozzle  shapes,  reducing  the  flow 
rate  through  the  propeller,  are  presented.  The  purpose  of  this 
reduction  is  to  prevent  the  occurrence  of  cavitation. 

The  investigation  was  carried  out  under  Contract  No.  N62558 
_3?C.O  of  the  U.S.  Department  of  the  Navy,  David  Taylor  Model  Basin, 
through  its  U.S.  Navy  European  Research  Contract  Program. 

-A  theory  is  described  for  the  numerical  calculation  of 

systematic  series  of  nozzles  shapes. | 

/ 

Three  nozzles  were  selected,  each  designed  for  operating 

at  the  same  thrust  coefficient  CT  ,  (  cT  =  o.gx  ) 

but  for  a  different  ratio  T  between  propeller  thrust  and  total 

thrust  (  T  =  1.04;  1.18;  1.36).>VThe  shape  of  the  nozzles  was 

chosen  to  produce  ram  pressures  at  the  propeller  plane  in 

order  to  delay  the  onset  of  cavitation.  Tunnel  experiments  were 

carried  out  with  the  above  nozzles. 

The  method  of  calculation  for  the  performance  characteristics 
of  nozzle  propeller  systems  gives  results  which  compare 
favourably  with  the  experimental  results. 
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1 .  Introduction. 

The  ducted  propeller,  invented  by  Kort  over  30  years  ago, 
is  now  extensively  used  in  cases  where  the  ship  screw  is 
heavily  loaded  (e.g.  tugs,  towboats  etc.).  Since  the  duct 
increases  the  flow  rate  through  the  propeller,  the  latter 
operates  at  a  more  favourable  loading. 

The  duct  itself  will  generally  produce  a  positive  thrust. 

The  application  of  this  kind  of  ducted  propellers  has  extensively 
been  dealt  with  in  the  literature.  The  theoretical  investigations 
of  Horn  [1]  ,  Dickmann  and  Weissinger [2]  ,  [3]  and  the 
systematic  experimental  investigations  of  van  Manen  [^1  ,  [5]  , 

may  be  mentioned  in  particular. 

The  range  of  applicability  of  the  ducted  propeller  may  be 
extended  since  the  duct  can  also  be  used  to  reduce  the  flow 
rate  through  the  propeller. 

This  second  type  of  flow  is  used  if  retardation  of  propeller 
cavitation  phenomena  is  desired.  In  this  case  the  duct 
reduces  the  flow  rate  through  the  propeller,  resulting  in  an 
increase  of  the  static  pressure  at  the  propeller  location. 

Ram  pressures  at  the  propeller  plane  are  obtained  if  the  mean 
static  pressure  at  the  propeller  plane  exceeds  the  static 
pressure  in  the  undisturbed  stream.  In  this  wa y  delay  of  the 
onset  of  propeller  cavitation  may  be  obtained.  The  duct 
itself  will  generally  produce  a  negative  thrust. 

Interest  has  recently  been  shown  in  the  application  of  the  second 
duct  type.  The  present  report  presents  the  results  of 
investigations  of  ducts,  suitable  for  building  up  "ram 
pressures"  at  the  propeller  location. 

The  investigation  covered  the  following  details: 

1.  Based  on  one-dimensional  momentum  considerat: ons,  expressions 
for  the  ideal  efficiency  of  the  propulsion  system  and  the  mean 
static  pressure  at  the  propeller  plane  were  derived. 

In  addition,  a  relation  between  total  thrust,  the  ratio 
between  propeller  thrust  and  total  thrust  and  the  angle 
.  between  nozzle  profile  and  shaft  line  is  given. 


2.  Representing  the  propeller  by  a  uniformly  loaded  actuator 
disk  rotating  with  finite  angular  velocity  and  the  nozzle 
by  distributions  of  ring  vortices,  sources  and  sinks,  a 
theory  was  developed  to  calculate  the  flow  field.  Starting 
from  the  flow  field,  expressions  for  the  propeller  torque 
and  thrust,  the  thrust  on  the  nozzle,  the  static  pressure 
along  the  nozzle  and  the  shape  of  the  camber  of  the  nozzle 
were  derived. 

3.  Systematic  calculations  were  carried  out  with  the  above 
theory.  The  relation  between  the  total  thrust  coefficient 
(KT  )  and  the  advance  coefficient  was  calculated  for  a 
number  of  nozzle  shapes.  In  addition  the  ratio  between 
propeller  thrust  and  total  thrust,  the  efficiency,  the 

mean  static  pressure  at  the  propeller  and  the  minimum  static 
pressure  at  the  exterior  surface  of  the  nozzle  -were  calculated. 

4.  Three  nozzles  were  selected,  each  designed  for  operating  at 
the  same  thrust  coefficient  CT,  but  for  a  different  ratio 
between  propeller  thrust  and  total  thrust.  The  shape  of  the 
nozzles  was  chosen  in  such  a  way,  that  "ram  pressures"  were 
expected  at  the  propeller  plane. 

Cavitation  tunnel  experiments  were  carried  out  with  tie  above 
nozzles  and  with  a  practical  nozzle  shape.  Five  screw 
models  have  been  tested  in  combination  with  each  of  the 
nozzles  .  The  propeller  torque,  the  propeller  thrust  and  the 
thrust  on  the  nozzle  were  measured. 

Flow  observations  were  performed.  Separation  phenomena  and 
cavitation  inception  at  the  surface  of  the  nozzle  were 
recorded. 

5.  The  validity  of  the  approximate  method  for  the  calculation 

of  the  nozzle-propeller  performance  was  tested,  by  a  comparison 
with  the  experimental  results. 


2.1.  Momentum  considerations. 

Insight  into  the  most  important  properties  of  ducted 
propellers  can  he  obtained  from  momentum  considerations. 
Figure  1  shows  diagrammatically  the  simplified  system  by 
which  the  ducted  propeller  can  be  replaced. 

The  K.y.z  body  axis  system  is  a  right  handed  orthogonal 
triad  with  its  origin  in  the  centre  of  the  propeller  plane. 
X  is  positive  in  the  direction  of  the  uniform  stream 
velocity  U. . 

The  propeller  was  replaced  by  a  uniformly  loaded  actuator 
disk.  The  tangentially  induced  velocities  were  neglected. 

The  nozzle  extends  from  x  =  _c/  to  *  =  with  a  radius  "R 

at  x  =  o  .  The  propeller  was  situated  in  the  middle  of 
the  nozzle. 

The  tobal  thrust  T  acting  on  the  fluid  owing  to  the 
working  propeller  and  no  the  nozzle  is: 

T  =  (u-t-u,)  (2.1  -  1) 

applying  the  momentum  theorem  over  the  control  volume 
given  in  Figure  1 . 

U,  and  Ux  are  the  mean  additional  velocities  in  the 
slip  stream  at  X=o  and  x-=od. 

A  denotes  the  actuator  disk  area  at  x=o. 

The  trust  Tp  developed  by  the  propeller  can  be  obtained 
by  calculating  the  pressure  difference  across  the 
actuator  disk  with  Bernoulli's  equation, 

"P.,!  +  j?(U  +  U.)2’  ijU* 

+  4KU+UJ 

and  thus, 


>(2.1-2) 


rp  =  A  =  <5  uu*A (i+  ^±) 


(2.1-3) 
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where  T;,  ■  and  T|A  denote  the  static  pressures 
upstream  and  downstream  of  the  actuator  disk,  respectively, 
and  denotes  the  static  pressure  for  upstream  and  for 
downstream  of  the  propeller. 

The  kinetic  energy  El  lost  in  the  propeller  slip  stream 
is  given  by, 

El  i y(u+u,) A (2.1-4) 

Fence  the  ideal  efficiency  yL  of  the  propulsion  device 
is  defined  by, 

^=-“1— = - 4=r  (2.1-5) 

UT+  £  i  +  VI+-VCt 

where 


A 

x  = 

T 

The  mean  value  of  the  static  pressure 
"p  can  be  calculated  from: 

"P  _  4-T,a 

meon  —  r 


>  (2.1-6) 

4  ' 

at  the  actuator  disk 


thus 


C  —  '^nnon  ~EST  _ T/^- _ 

Pmcon~  "  *  l_l+/rnEc?. 


(2.1-7) 

The  efficiency  yL  and  the  mean  static  pressure  coefficient 
at  the  actuator  disk  ce  are  plotted  as  a  function  of 
the  ratio  T  between  propeller  thrust  and  total  thrust, 
wil'h  the  total  thrust  coefficient  cT  as  parameter,  in 
the  Figures  2  and  5  respectively. 


V 


1?. 
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It  follows  from  these  Figures  that  the  efficiency  of 
the  propulsion  device  decreases  and  the  mean  value  of  the 
static  pressure  at  the  propeller  plane  increases  with  in¬ 
creasing  value  of  the  ratio,?. Ram  pressures  at  the  propeller 
plane  are  only  built  up  when  X  exceeds  1.0.  Consequently 
a  negative  thrust  is  acting  on  the  duct  in  that  case. 

The  total  thrust  coefficient  cT  and  the  ratio  T  between 
propeller  thrust  and  total  thrust  can  approximately  be 
calculated  as  a  function  of  the  angle  of  incidence  of  the 
nozzle  profile  in  the  way  as  described  by  Chen  [7]  »  [8] 

The  radial  velocity  UA  ,  induced  by  the  propeller  at  the 
nozzle,  is  principally  a  function  of  the  actuator  disk 
loading.  The  results  given  by  lim  and  Chen  [9]  < 
enable  the  calculation  of  the  mean  value  of  the  radially 
induced  velocity  along  the  duct  as  a  function  of  the 
propeller  thrust  coefficient  CTp 
between  the  nozzle  chord  length 
radius  a  . 


and  the  r«v,io 
c  and  the  propeller 


(2.1-8) 


The  propeller  may  be  considered  as  operating  in  open 
water  with  an  equivalent  uniform  stream  velocity 
u+un=  u  +  u,-^  .  Therefore  equation  (2.1-8) 

becomes: 


— ^  (  ctp  ;  J|  ) 
u+u_  *  p  -R 


(2.1-9) 


where  un  denotes  the  mean  axial  velocity  induced  by 
the  nozzle  at  the  propeller  location. 

The  sectional  lift  of  the  nozzle  is  according  to  the 
two  dimensional  theory  equal  to  (see  for  sign  conventions 
Figure  4) : 
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^  cCUo< 


(2.1-10) 


where  CL<<  and  <*  denote  the  sectional  lift  coefficient 
and  the  angle  of  attack  of  the  nozzle  profile  respectively 
In  the  case  of  symmetrically  loaded  ducts  there  are  no 
radial  velocities  induced  hy  the  nozzle  in  the  mid  plane 
of  the  nozzle.  Consequently  the  thrust  coefficient 
of  the  nozzle  cT  can  he  readily  obtained  from 
equations  (2.1-9)  and  (2.1-10). 


_  X 


ijuNrs* 


=  a-l  Cl-<*  ^  (CV-|) 


(2.1-11) 


Furthermore , 


CT  =  (*--c)CT 


(2.1-12) 


From  equations  (2.1-9),  (2.1-11)  and  (2.1-12)  it 


follows  that 


0-tQct _ 

sir,[«RC^ 


v/here  ©  denotes  the  angle  between  the  nozzle  p.ofile  and  th 
shaft  line.  The  relation  between  the  thrust  coefficient 
Ct  and  the  ratio  t  =  is  plotted  in  Figure  5  for  some 

T* 

values  of  G  and  for  the  case  £  -  i.a  and 

— R 

Figure  5  shows  that  for  a  given  nozzle  shape,  the  thrust 
on  the  nozzle  increases  with  increasinf.  loading  of  the 
nozzle-propeller  system. 


itjiiii.-  i^pumt  yi. 
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2.2.  Representation  of  ducted  propellers  by  vortex 
distributions. 

The .calculations  on  the  ducted  propellers  are  phased  on 
the  following  assumptions.  The  ducted  propeller  moves 
steadily  forward.  The  forward  velocity  was  assumed  to 
be  sufficiently  large,  the  nozzle  loading  and  the  propeller 
loading  sufficiently  low  to  permit  the  application  of 
linearized  theory. 

The  calculations  on  ducted  propellers  may  be  classified  into 
two  types: 

(1)  the  direct  problem,  in  which  the  propeller  blade  form 
and  the  shape  of  the  duct  are  described, 

(2)  the  inverse  problem,  in  which  a  certain  combination 
of  blade  forces  and  duct  forces  are  given. 

Notable  v/ith  respect  to  the  direct  problem  is  the 
theoretical  tr ?atment  of  Morgan  [5]  . 

The  present  investigation  is  concerned  with  the 
inverse  problem. 

The  considered  ducted  system  consists  of  an  annular 
airfoil  of  finite  length  with  an  impeller  having  an 
infinite  number  of  blades.  The  mathematical  model  of 
the  geometrical  configuration  of  the  ducted  propeller 
can  be  represented  by  vortex  and  source  distributions 
as  summarized  in  Figure  6. 

The  propeller  is  considered  as  an  actuator  disk  which 
is  set  normal  to  the  free  stream.  It  is  driven  to 
rotate  in  its  own  plane  at  an  angular  velocity  go. 

The  propeller  flow  field  consists  of  helical  trailing 
vortices  starting  from  the  propeller  disk  at  hub  and 
nozzle  diameter. 

The  strength  of  the  trailing  vortices  is  a  function  of 
the  loading  of  the  propeller  disk  y  and  the  advance 
coefficient  3 . 
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Each  helical  trailing  vortex  line  lies  on  a  cylinder 
of  constant  diameter  (equal  to  huh  or  nozzle  diameter) 
and  "has  a  constant  pitch. 

The  radial  component  of  the  induced  velocity  exhibits  a 
logarithmic  singularity  at  the  periphery  of  1  he  propeller 
disk.  Since  logarithmic  singularities  are  integrable, 
the  stream  line  near  the  propeller  tip  is  continuous  hut 
with  infinite  slope. 

The  flow  around  the  nozzle  is  represented  by  a  distribution 
of  ring  sources  and  a  distribution  of  ring  vortices  along 
a  C7/Iinder  of  constant  diameter.  The  annular  airfoil  is 
axisymmetrical,  so  that  the  nozzle  has  no  trailing  vortices. 
The  nozzle  is  thus  replaced  by: 

(1)  A  bound  ring  vortex  distribution  v/ith  a  strength  equal 
to  zero  at  the  leading  edge  of  the  nozzle  and  equal 

to  the  stren.  th  of  the  circumferential  component 
of  the  helical  trailing  vortices  3t  the  propeller  disk. 
The  induced  radial  velocity  due  to  this  vortex 
distribution  has  a  logarithmic  singularity  at  the  _ 
propeller  plane  at  the  tip  diameter.  The  ring  vortex 
distribution  along  the  duct  has  been  chosen  in  such 
a  way  that  the  logarithmic  singularity  of  the  radial 
velocity  induced  by  this  vortex  distribution  and  the 
actuator  disk  compensate  each  other.  Consequently  the 
total  radialy  induced  velocity  has  a  smooth 
behaviour.  This  fact  simplifies  the  numerical 
calculations* 

(2)  A  source  and  sink  distribution  representing  the 
thickness  effect  of  the  nozzle. 

(3)  Continuous  bound  ring  vortex  distributions  v/ith  zero 
strength  at  the  leading  and  trailing  edges  of  the 
nozzle  and  v/ith  sinusoidal  shapes.  Other  vortex 
distributions  along  the. nozzle  can  be  built  up  by 
Fourier  synthesis. 

The  nozzle  will  have  shock  fT'ee  entry  because  the  total 
ring  vortex  stro:  gth  at  the  leading  edge  equals  zero. 
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The  range  of  angles  of  attack  within  which  this 
shock  -  free  entry  can  he  realized  if  the  loading  of  a 
given  ducted  propeller  system  is  varied  largely  depends 
on  the  shape  of  the  leading  edge  of  the  nozzle  profile. 
Shock-free  entry  occurs  on  thick  profiles  having 
rounded  leading  edges,  over  a  largo? range  of  angles 
of  attack  than  on  thin  profiles  having  sharp  leading 
edges. 

The  resulting  mathematical  model  is  summarized  in 
Figure  6,  where  wa  ;  wT  and  are  the  dimensionless 
velocities  (in  axial-,  circumferential-  and  radial 
direction,  respectively)  induced  hy  all  the  vortices 
and  sources.  The  vortex  strength  per  unit  area  of  the 
actuator  disk  is  denoted  hy  y  .  The  amplitudes  of  the 
sinusoidal  ring  vortex  distributions  along  the  nozzle 
are  denoted  hy  .  T  .e  source  and  sink  distribution 
along  the  nozzle  is  given  hy 

The  total  induced  velocities  can  he  calculated  according 
to  the  law  of  Biot  -  Savart.- 

The  propeller  flow  field  is  known  now  and  the  propeller 

thrust,  the  thrust  of  the  nozzle,  the  propeller  torque 

% 

and  the  efficiency  can  he  calculated  for  the  chosen  data 
of  design  parameters.  At  the  same  time  the  shape  of  the 
camber  line  of  the  nozzle  and  the  static  pressure  along 
the  nozzle  can  he  calculated. 

The  details  of  the  theory  are  not  dven  here. 

A  numerical  program  for  the  \.C.I.:.3.  -  digital  computer 
based  on  this  theory  was  made. 

Computations  were  carried  out  to  select  systematic 
series  of  nozzle  shapes. 
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3*  Calculation  of  systematic  series  of  duct  shapes. 

The  basic-design  parameters  of  ducted  propellers  are  (see 
Figure  ?): 

(1)  the  ratio  between  nozzle  length  and  propeller  diameter,  9. 

(2)  the  ratio  between  nozzle  length  in  front  of  propeller 
disc  and  tots'1  nozzle  length,  £. 

(3 )  the  ratio  between  hub  and  propeller  tip  diameter,  jL 

(4)  the  thickness  distribution  and  the  maximum  thickness  of 
the  nozzle  profile. 

(5)  the  advance  coefficient  3  . 

(6)  the  loading  of  the  propeller. 

(7)  the  loading  of  the  nozzle. 

Based  on  the  theoretical  analysis  described  in  the  previous 
section  calculations  were  carried  out  to  come  to  series  of 
nozzle  shapes. The  data  used  for  the  design  parameters  are: 

S.  3  0.6 
■y> 

2  =  0.5 
c. 

d  n  o.x 

3 

The  nozzle  profile  has  a  NASA  oois 
basic  thickness  form. 

The  loading  of  the  propeller,  the  loading  of  the  nozzle  and 
the  advance  coefficient  were  systematically  varied. 

The  numerical  results  are  presented  in  Figure  2  and 
graphically  represented  in  the  Figures  8  through  12. 


The  effect  of  the  ratio  T  between  propeller  thrust  and  total 
thrust  on  efficiency  for  some  values  of  the  thrust  coefficient 
cT  and  for  is  shown  in  Figure  2. 

The  ducted  propeller  corresponds  for  the  case  h=o  with  an 
actuator  disk  enclosed  by  a  nozzle  and  rotating  at  an 
infinite  angular  velocity. 

The  results  of  the  momentum  considerations  described  in 
section  2.1  are  presented  in  Figure  2. 

The  efficiencies  calculated  by  these  two  methods  agree 
very  well  for  small  values  of  r  .  The  differences  between 
the  efficiencies  for  large  values  of  T  can  be  explained  by 
the  thickness  effect  of  the  nozzle  which  is  only  taken  into 
account  in  the  vortex  theory.  The  effect  of  the  nozzle  cn  the 
flow  at  the  propeller  is  in  that  theory  taken  into  account 
in  a  more  thorough  way. 

The  differences  between  the  curves  in  Figure  2  show  that  for 
large  ratios  between  pro. -oiler  thrust  and  total  thrust, 
thin  nozzle  profiles  are  more  suitable  with  respect  to 
efficiency  «• 

The  efficiency  of  the  propulsion  system  always  diminished  with 
increasing  ratio  t  .  However,  it  is  noted  that  for  lightly  loaded 
systems,  the  decrease  in  efficiency  is  small. 

The  shape  of  the  camber  line  of  a  nozzle  -SC*)  is  completely 
determined  by  the  strength  of  the  vortices  along  the  nozzle 
y,  the  ratio  between  the  vortex  strength  at  the  propeller 
plane  ^  ana  the  advance  coefficient  b  and  by  the 
geometry  of  the  system  (2;£  ;d  ;  thickness  distribution 
and  maximum  thickness  of  the  nozzle). 

sex')  =  f  C  y.  >  5  >  qeoMfiTRy). 

The  propeller  is  represented  by  an  actuator  disk  rotating  at 
an  infinite  angular  velocity  if  the  undisturbed  stream  velocity 
u.  is  assumed  to  be  constant  and  the  advance  coefficient  b 
becomes  zero. 
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The  ratio  l  is  kept  constant,  thus  the  vortex  strength  at  the 
propeller  goes  in  the  same  way  to  zero  as  the  advance 
coefficient  3. 

The  total  thrust  coefficient  cT  and  the  ratio  between  propeller 
thrust  and  total  thrust  T  are  denoted  by  cTo  and 
if  3=o  .  Thus  the  shape  of  a  nozzle  is  also  completely 
determined  by  cT<  ;  ~co  and  the  geometry  of  the  ducted 
propeller  system. 


SCx)  ■=  f  (c-To  ;  ;  qeowETRy) 

Calculations  of  the  efficiency  ,  the  total  thrust 
coefficient  KT  and  the  ratio  between  propeller  thrust  and 
total  thrust  at  various  advance  coefficients  3*, 
we re  made  for  a  number  of  nozzle  shapes  determined  by 
cT  and  t  . 

In  addition,  the  mean  value  of  the  static  pressure  at  the 
propeller  plane  and  the  minimum  static  pressure  at 

che  exterior  surface  of  the  nozzles  T  •  were  calculated 

“in 

for  the  nozzles  considered. 

The  following  non-dimensional  pressure  coefficients  are 
introduced. 


_  mean 


-'Ro 


P  min 


where  "•%  and  are  the  static  pressure  and  the 

dynamic  pressure  of  the  undisturbed  stream  respectively. 


The  numerical  results  are  presented  in  the  Figures  8 
through  12.  The  nozzle  shapes  considered  are  tabulated  below. 
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It  appears  that  the  efficiency  17.  of  a  ducted  propeller 
decreases  for  increasing  advance  coefficient  iJ. 

This  phenomenon  can  be  explained  by  the  losses  due  to 
rotation.  The  kinetic  energy  in  the  propeller  slip  stream, 
which  is  lost,  may  be  split  up  into  the  losses  due  to  the 
axial  acceleration  of  the  fluid  and  those  due  to  the 
rotation  of  the  fluid.  The  losses  lue  to  the  axial 
acceleration  are  independent  of  the  advance  coefficient  iJ , 
those  due  to  the  rotation  of  the  fluid  are  equal  to  zero 
for  3=o  and  increase  with  increasing  b  . 

Application  of  pre -turning  vanes  becomes  important  in  the 
case  of  heavily  loaded  systems  if  the  gain  in  rotational 
efficiency  is  larger  than  the  efficiency  loss  due  to 
friction. 

The  thrust  coefficient  cT  of  the  ducted  propeller  always 
diminishes  with  increasing  advance  coefficient  3  .  This 
phenomenon  can  also  be  explained  by  the  rotation  of  the 
fluid. 

The  minimum  static  pressure  at  the  exterior  surface  of  the 
nozzle  is  almost  independent  of  the  total  thrust  coefficient 
cT  .  V.'ith  increasing  ratio  x  between  propeller  thrust  and 
total  thrust,  the  minimum  static  pressure  decreases. 
Consequently  the  risk  of  cavitation  at  the  exterior  surface 
of  the  nozzle  increases  in  that  case. 


% V'  '  - 

"<  ■'  - 

' 
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The  shape  of  the  nozzles  for  the  experiments  were  selected  in 
such  a  way  that  "ram  pressures"  are  expected  at  the  propeller 
plane.  Thus  the  ratio  should  be  above  1.0.  The  nozzles 

are  determined  for  relatively  high  ship  speeds. 


The  data  used  for  the  nozzles  are: 


nozzle 


O.gX 

Xn  =  i  .oA 

O.g?. 

i.ift 

°sx 

i.s6 

The  nozzle  shapes  are  presented  in  Figure  13  and  tabulated  in 
table  1 . 


The  ideal  efficiency  ,  the  thrust  coefficient  KT  > 

and  the  ratio  x  -Tp  at  various  advance  coefficients  'd  are 

*r 

presented  in  Figure  14.  In  addition  the  pressure  coefficients 


and 


mean 


P  mm 


are  given. 


The  experiments  were  carried  out  in  the  N.S.M.B.  cavitation 
tunnel  no.  1  having  a  90  cm  x  90  cm  closed  test  section  and 
a  uniform  flow.  Five  different  screw  models  were  tested 
in  combination  with  the  nozzles  (j§)  (26)  (30)  (3j)  and  (5$) 

The  results  of  earlier  investigations  by  Van  Manen  into 
screws  in  nozzles  are  given  in  the  N.S.fil.B.  publications  (V] 

[5] 

Nozzle  No.  (b)  is  recommended  by  the  N.S.M.B.  for  practical 
purposes.  This  nozzle  increases  the  flow  rate  through  the 
propeller  disk. 

Later  on  it  was  tried  to  build  up  a  ram  pressure  inside  the 
nozzle  by  a  suitable  choise  of  the  camber  of  the  nozzle  profile.* 
A  practical  shape  of  a  nozzle  was  designed,  indicated  by 
No .  (26)  .  Nozzle  No.  (»6)  has  a  cylindrical  inner  wall  so  that 

the  screw  can  arbitrarily  be  located  in  the  nozzle  with  a 
constant  clearance  between  blade  tip  and  nozzle. 

Nozzle  No.  'so)  (5?)  and  (32)  are  calculated  with  the  aid  of 
the  theoretical  analysis  described  in  the  previous  sections. 

The  mean  lines  and  the  thickness  distributions  of  all  the 
mentioned  nozzle  shapes  are  presented  in  Figure  13. and 
tabulated  in  table  1 . 

The  experiments  were  all  carried  out  v/ith  a  series  of  five 
b'laded  Kaplan  type  screws.  The  most  important  characteristics 
ire  presented  in  the  table  below.  Further  particulars  of  the 
screw  models  are  given  in  the  Figures  15  through  19* 
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table 

Diameter 

3 

240  mm 

Number  of  blades 

z 

5 

Pitch  ratio 

p/D 

1.0-1 .2-1 .4-1. 6-1. 8 

Pitch  distribution 

uniform  near  the  tip, 
decreasing  to  the  hub. 

Blade  area  ratio 

B.A.R. 

100  % 

Blade  outline 

Kaplan  type 

Blade  section 

NASA  16  -  parabolic 

camber  line  - 

flat  face. 

Total  number  of  screws 

5 

Propellers  indicated  by 

3527,  5528,  5529, 

5550,  5551. 

The  screws  were  designed  in  combination  with  nozzle  No.  ^9) 

The  pitch  distribution  depends  on  the  nozzle  induced  velocities 
at  the  location  of  the  screw  and  on  the  racial  load  distribution 
on  the  screw. 

>  The  screws  were  located  in  the  nozzles  with  a  uniform  tip 

'•  clearance  of  1  mm. 

Cavitation  tunnel  experiments  were  carried  out  at  constant 
rpm  of  the  screw  models  ,  at  speeds  of  advance  covering  a  slip 
range  from  4o  -  100  %  and  at  various  cavitation  numbers.  The 
torque  and  thrust  on  the  screw  and  the  thrust  on  the  nozzle 
were  recorded  by  means  of  strain  gauges. 

The  total  thrust  coefficient  KT  ,  the  torque  coefficient 
the  efficiency  ^  and  the  ratio  t  between  propellqr  thrust  and 
total  thrust  were  all  calculated  as  functions  of  the  advance 
coefficient  3 .  For  the  nozzles  (ji)  ,  (2&)  ,  (0)  ,  (sB)  and 

the  experimental  results  are  presented  in  the  Figures '^0 
through  24  respectively. 


In  addition,  the  mean  static  pressure  coefficient  CPme/an 
was  derived  from  the  measured  thrust  coefficient  cT  and  the 
ratio  *t  with  the  aid  of  Figure  3*  The  result  is  also  presented 
in  the  Figures  20  through  24* 

The  relation  between  the  thrust  coefficient  cT  and  the 
ratio  x  for  the  various  nozzle  shapes  is  given  in 
Figure  5»- 

The  results  of  the  experiments  are  synthesized  in  a  3p_8 
diagram,  (see  Figure  25).  The  parameters  3  and  6  are 
defined  by: 

k  * 

=  33.08  TQ 
p  3%. 

S  -  lo>-*» 
b 

Cavitation  observations  were  made  during  all  the  experiments. 
Figure  26  shows  the  inception  lines  for  bubble  cavitation  at 
the  exterior  surface  of  the  nozzles. 

The  cavitation  number  a~0  is  plotted  to  a  base  of  total 
thrust  coefficient  cT  .  The  cavitation  number  cr  is' 
defined  by 


where  denotes  the  vapour  pressure  of  the  water. 

The  pressure  coefficient  cpmin  ,  determining  the  minimum 
static  pressure  ai  the  exterior  surface  of  the  nozzle,  was 
calculated  for  the  various  nozzles. 

The  numerical  results  are  presented  in  Figure  26. 

It  is  noted  that  the  results  of  the  thrust,  torque  and  velocity 
measurements  were  not  corrected  for  the  wall  effects  of  the 
cavitation  tunnel. 
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5.  Discussion  of  theoretical  and  experimental  results. 

Certain  simplifying  assumptions  were  made  in  the  development 
of  the  theo.-.y  for  the  numerical  calculations  of  nozzle  shapes. 
/  The  validity  of  the  method  must  he  tested  by  a  comparison 
with  measurements. 

The  experimentally  obtained  relation  between  the  thrust 
coefficient  KT  and  the  advance  coefficient  "3  for  the  nozzles 
No.  (0)  ,  (0)  and  in  combination  with  the  various 

screws,  is  given  in  the  Figures  22  through  24  respectively. 

The  theoretically  obtained  relations  between  KT  and  T 
for  which  the  nozzles  No.  ,  (^)  and  were  designed 

are  presented  in  Figure  14.  A  nozzle  screw  combination  meets 
the  requirements  of  tne  theory  if  the  theoretical  design  curve 
of  the  nozzle  and  the  measured  thrust  curve  of  the  nozzle  - 
screw  combination  intersect  at  the  dcsrgn  b  . 

The  point  of  intersection  of  the  various  nozzle  -  screw 
systems  are  given  in  the  following  table. 

The  measured  torque  coefficient  ,  the  efficiency  7  and 
the  ratio  t  are  tabulated  aswell. 


Nozzle  No. 

screw 

% 

Kt 

d 

Wq 

7 

X 

3527 

1.0 

0.193 

0.75 

0.0379 

0.60 

1.05 

3528 

1.2 

0.262 

0.88 

0.0588 

0.63 

1.05 

3529 

1.4 

0.337 

1.00 

0.0834 

0.65 

1.05 

3530 

1.6 

0.409 

1.11 

0.1101 

0.67 

1.03 

3531 

1.8 

0.490 

1.24 

0.1398 

0.69 

1.01 

© 

3527 

1.0 

0.209 

0.79 

0.0430 

6. 61 

1.19 

3528 

1.2 

0.279 

0.92 

0.^672 

0.61 

V" 

C\J 

• 

3529 

1.4 

0.342 

1.03 

0.0941 

0.60 

1.22 

3530 

1.6 

0.404 

1.13 

C.1244 

0.59 

1.23 

3531 

1.8 

0.464 

1.23 

0.1596 

0.58 

1.25 

© 

3527 

1.0 

0.223 

0.82 

0.0558 

0.53 

1.17 

3528 

1.2 

0.300 

0.97 

0.0823 

0.56 

1.18 

3529 

1.4 

0.360 

1.10 

0.1126 

0.57 

1.18 

3530 

1.6 

0.420 

1.22 

0.1462 

0.55 

1.24 

3531 

1.8 

0.456 

1.30 

0.1836 

0.51 

1.20 

The  above  results  afford  the  possibility  of  comparison  with  the 
predicted  performance . 

The  overall  efficiency  of  the  ducted  propeller  may  be  expressed 
as  a  product  of  the  ideal  efficiency  y.  and  the  efficiency 
factor  due  to  profile  drag  and  frictional  drag  of  the 
propeller  blades  and  the  nozzle. 
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For  low  advance  coefficients  ^  ,  the  screw  is  represented  in 
the  theory  by  ai  actuator  disk  rotating  at  a  high  angular 
velocity.  Consequently  the  efficiency  loss  of  the  screw  due  to 
friction  becomes  large  and  the  efficiency  factor  fj.  tends 
to  zero.  The  frictional  losses  decrease  with  increasing  hi. 

Figure  14  shows  that  the  agreement  between  the  calculated  values 
of  y  and  those  obtained  by  the  tests  iS  good  if  the  above 
considerations  with  respect  to  the  frictional  losses  are 
taken  into  account. 

It  is  also  seen  from  Figure  14  that  the  calculated  values  of  the 
ratio  T  between  propeller  thrust  and  total  thrust  agree  very 
well  with  the  test  results. 

Flow  observations  were  performed  during  the  cavitation  tunnel 
experiments.  Separation  phenomena  at  the  exterior  or  the 
interior  surface  of  the  nozzles  No.  ($d)  ,  (|l)  and  (££) 

we re  not  observed  over  the  range  of  considered  water  speeds. 

The  shapes  of  the  cambers  of  nozzles  No.  and  No.  (0)  are 

geometrically  almost  equal,  (see  Figure  13)*  The  cambers  are 
only  different  at  the  leading  edges  of  the  nozzle  profiles. 

Flow  separation  was  observed  at  the  interior  surface  of  nozzle 
No.  (^6)  during  the  cavitation  tunnel  tests,  This  phenomenon 
can  be  explained  by  the  unfavourable  shape  of  the  nozzle  at 
the  leading  edge. 

The  slope  of  a  streamline  passing  through  the  propeller  tip 
of  a  propeller  in  open  water  is  directly  proportional  to  the 
propeller  thrust  coefficient  cTp  •  The  thrust  coefficient  cTh 
is  primarily  determined  by  the  thrust  coefficient  of  the 
propeller  cTp  and  the  nozzle  geometry. 

Hence  it  is  expected  that  there  exists  a  fixed  relation 
between  the  thrust  coefficient;  c-  and  the  ratio 
of  a  nozzle  which  relation  is  independent  of  the  advance 
coefficient  3  or  of  the  screw  considered. 
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The  experimentally  obtained,  relations  between  CT  and  x  for  the 
nozzles  No.  ,  (^)  ,  ,  (0)  and  (^2)  ,  in 

combination  with  the  various  screws  are  given  in 
Figure  5  .  The  results  confirm  the  above  supposition. 

The  thrust  coefficient  CT  and  the  ratio  x  between  propeller 
thrust  and  total  thrust  are  approximately  calculated  as  a 
function  of  the  angle  ©  between  the  nozzle  profile  and  the 
propeller  shaft  line  in  section  2.1.  The  numerical  results  are 
also  presented  in  Figure  5  » 

The  agreement  between  the  calculated  relations  of  cT  and  x  =  1e 
and  those  obtained  by  the  tests  is  good  in  the  neighbourhood 
of  the  design  points  of  the  various  nozzles. 

Inception  lines  for  bubble  cavitation  at  the  exterior  surface 
of  the  nozzles  are  presented  in  Figure  26. 

A  comparison  between  the  inception  lines  of  the  various  nozzles 
shows  that  for  equal  thrust  coefficients  CT  ,  the  inception  of 
cavitation  shifts  co  higher  at  increasing  negative  loading  on 
the  duct  or  at  correspondingly  increasing  ratio  . 

The  numerically  obtained  maximum  cavitation  numbers 
are  also  presented  in  Figure  26  for  the  various  nozzles. 

The  calculated  cavitation  numbers  show  the  same  tendency  as  the 
measured  inception  lines. 

The  differences  between  the  calculated  and  the  measured 
cavitation  numbers  may  be  explained  by  the  lack  of  exact  data 
on  cavitation  inception. 

Thus  far,  no  remarks  were  made  on  the  cavitation  characteristics 
of  the  screw  series  considered.  The  screws  were  designed  in 
combination  with  nozzle  (1^)  .  The  pitch  distribution  depends 

on  the  nozzle  induced  velocities  at  the  propeller  plane  and  on  the 
radial  load  distribution  of  the  screw.  The  distribution  of  the  nozzle 
induced  velocities  at  the  propeller  plane  is  rather  different  for 
the  nozzles  (^)  ,  (f^)  ,  (||)  ,  and  .  Consequently  the  screws  have  bad 
characteristics  with  respect  to  cavitation  inception  if  these  are 
t'ested  in  combination  with  the  nozzles  ,  (£l)  and  (^2)  .  Therefore 
no  cavitation  observations  of  the  screws  are  reported. 


- v 
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6.  G one lus ions. 

1.  The  ideal  efficiency  of  a  ducted  propeller  system 
decreases  with  increasing  ratio  T  between  propeller  thrust 
and  total  thrust.  However,  it  is  noted  that  for  lightly 
loaded  systems,  the  decrease  in  efficiency  is  small. 

2.  Especially  in  the  case  of  decelerating  nozzles,  it  is 
recommendable  to  choose  thin  nozzle  profiles  due  to  the 
unfavourable  effect  of  thickness  of  the  nozzle  on 
efficiency. 

J.  The  mean  static  pressure  at  the  propeller  plane  increases 
with  increasing  ratio  X  betv/een  propeller  thrust  and  total 

,  thrust.  "Ram  pressures"  at  the  propeller  plane  are  only 
built  up  if  x  becomes  larger  than  1.0,  consequently  a 
negative  thrust  is  produced  by  the  duct  in  that  case. 

4.  The  method  of  calculation  for  the  performance  characteristics 
of  a  nozzle  propeller  system  developed  in  this  study  gives 
results  which  compare  favourably  with  the  experimental 
results.  It  is  recommended  for  a  possible  extension  of  the 
research  to  determine  in  the  f'.rst  instance  the  effects  of 
the  various  design  parameters  ( 4  ;  £  and  so  on)  on  the 
performance  characteristics,  with  the  aid  of  the  theory. 

5.  The  application  of  pre-turning  vanes  becomes  important  if  the 
gain  in  rotational  efficiency  is  larger  than  the  efficiency 
loss  due  to  friction  of  the  vanes.  This  will  be  the  case  if 
the  ducted  propeller  system  is  heavily  loaded  and  the 
angular  velocity  of  the  propeller  is  relatively  small. 

The  losses  due  to  the  rotation  of  the  fluid  can  be  isolated 
with  the  aid  of  the  theory.  Hence,  the  losses  due  to 
rotation  and  the  loss  due  to  the  friction  of  pre-turning  vanes 
may  be  easily  compared. 

6.  The  static  pressure  at  the  exterior  surface  of  the  nozzle 
decreases  if  the  static  pressure  at  the  propeller  plane 
increases.  However,  the  reduction  in  static  pressure  at  the 
exterior  nozzle  surface  is  small  in  comparison  with  the 
gain  in  pressure  at  the  propeller  plane. 
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7.  It  follows  from  theory  and  experiment  that  there  exists 
a  fixed  relation  between  the  thrust  coefficient  cT  and  the 
ratio  t =  "5  of  a  nozzie,  which  relation  depends  neither 
on  the  advance  coefficient  3  nor  on  the  screw  considered. 
This  fact  gives  in  an  easy  way  information  on  the  range 
of  application  of  the  nozzle. 
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BLZ. 

NC. 

_ _ _ _ 

X 

nozzle  30 

lozzle  31 

Sb% 

nozzle  3- 

sm/l 

0 

+ 

0 

• 

0 

_A 

ro 

-  0.0073 

-  0.0232 

0.083 

+  0.0130 

+  0.0008 

-  0.0096 

0.16? 

+  0.0139 

+  0.0067 

+  0.0006 

0.250 

+  0.0128 

+  0.0093 

+  0.0067 

0.333 

+  0.0102 

+  0.0091 

+  0.0085 

0.41? 

+  0.0059 

+  0.0060 

+  0.0063 

0.500 

+  0 

+  0 

+  0 

0.583 

-  0.0081 

-  0.0096 

-  0.0112 

0.667 

-  0.0186 

-  0.0229 

-  0.0276 

0.750 

-  0.0316 

-  0.0404 

-  0.0496 

0.833 

-  0.04-76 

-  0.0622 

-  0.0773 

0.017 

-  0.0661 

-  0.0881 

-  0.1106 

1.000 

-  0.0874 

-  0.1182 

-  0.1498 

Table  1 
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Figure  3.  Mean  static  pressure  at  propeller  plane  of 
propeller  nozzle  system. 
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1.  Screw  disk  with  bound  vortices. 

2.  Helical  trailing  vortices. 

3.  Discontinuous  ring  vortex  distribution  along  nozzle 

4.  Continuous  ring  vortex  distribution  along  nozzle 

[Y(x)=Ymsinm©]. 

5.  Source  and  sink  distribution  along  nozzle. 


Figure  6.  Mathematical  model  of  ducted  propeller  with 
infinite  number  of  blades. 
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Figure  1G.  Characteristics  ol‘  a  systematic  series  of  nozzle  shapes. 
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Figure  11.  Characteristics  of  a  systematic  series  of  noxale  shapes. 
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Figure  12.  Characteristics  of  a  systematic  series  of  nozzle  shapes 
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Figure  t4.  Characteristics  of  nozzles  selected  for  experiments. 
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Fjg. 15,16, 17,18 and  19:  Particulars  of  screw  models  3527 
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id  IS.  Particulars  of  screw  models  3527  through  3531. 


Figure  20.  Experimentally  obtained  characteristics  of  nozzle  NS 
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Figure  21.  Experimentally  obtained  characteristics  of  nozzle 


Figure  22.  Experimentally  obtained  characteristics  of  nozzle 
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Figure  23.  Experimentally  obtained  characteristics  of  nozzle 
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Figure  26.  Inception  lines  for  bubble  cavitation  at  the  exterior  surface  of  the  nozzles 


